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Abstract. We derive explicit formulas for the normal ordering of powers of arbitrary 
monomials of boson operators. These formulas lead to generalisations of conventional Bell 
^ | , and Stirling numbers and to appropriate generalisations of the Dobinski relations. These 

new combinatorial numbers are shown to be coherent state matrix elements of powers of 
the monomials in question. It is further demonstrated that such Bell-type numbers, when 
considered as power moments, give rise to positive measures on the positive half-axis, which 
in many cases can be written in terms of known functions. 

OO ! 

q ■ The Standard boson commutation relation [a, af] = 1 can be realised by identifying 

\ formally a = j- and al = x, since [4-, x] = 1. In the present note we shall use both forms. 
Integer sequences arise naturally when considering the action of (x d/dx) n on f(x) which in 
q ■ general can be written as 

ÇL; (x ^Tf(x) = £ S ( n > k)x k (d/dx) k f(x) (1) 

I fc — 1 

or, alternatively [[[]] 

Cf (ala) n = Y,S(n,k)(a^) k a k . (2) 

> ! k=l 

The Stirling numbers ofthe second kind S(n, k) appearing in Eqs.([l]) and (0) have been known 
for over 250 years 0. Eq.(@) exemplifies the normal ordering problem, that is, finding the 
form of (a)a) n with the powers of a on the right. Although explicit expressions for S(n, k) 
are known [Q], of particular interest here are the Bell numbers B{n) given by the sums 

n 

B(n) = Y,S(n,k), n = l,2,... (3) 

k=l 

with 5(0) = 1 by convention (S(n, 0) = <5 ní0 ). A closed-form expression for B(n) can be 
found by considering the action of (xd/dx) n on f(x) = e x giving: 

oo r,n n 

(V^) H Tf = £ S ( n > k ) xk i W 
k=o K - k=l 

which for x = 1 reduces to 

oo i n n 

(l/e)Y,T^ = T,S(n,k) = B(n). (5) 

k=0 K - k=l 

Equations (gj) and (Q) are the celebrated Dobinski formulas ([0],[0],[0],[0]) which have been 
the subject of much combinatorial interest. Eq.@ represents the integer B(n) as an infinite 
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series, which is however not a power series in n. An immediate consequence of Eqs.(Q) 
and @ is that B(n) is the n-th moment of a (singular) probability distribution, consisting of 
weighted Dirac delta functions located at the positive integers (the so-called Dirac comb) : 



B(n) = / x n W(x)dx, n = 0, 1, . . . (6) 
J o 

where 

W(x) = (l/e)± 6 -^^. (7) 

The discrete measure W(x) serves as a weight function for a family of orthogonal polynomials 
C^(x), the Charlier polynomials [|p. The exponential generating function (EGF) of the 
sequence B{n) can be obtained from Eq.@ as 



e 



oo \n 
( e *-l) v^„, >A 



k=0 



This equation is related via Eq.(^) to a formula giving the normal ordered form of e Aíïta |]7p,[H] 

The symbol jV denotes normal ordering, while :: refers to normal-ordering with neglect of 
non-commutativity. We stress that in the derivation of Eq.@ in [0] and ^ no use has been 
made of the Stirling and Bell numbers. It may readily be seen that Eq.(§) is the expectation 
value of Eq.(0) in the coherent state \z) defined by a\z) = z\z) at the value \z\ = 1, using 
Eqs.(@) and @. This circumstance has been used recently to re-establish the link between the 
matrix element (z|e Aata |z) and the properties of Stirling and Bell numbers [§]. 

The purpose of this note is to show that the above results on functions of a)a may be 
extended to functions of (a^) r a s , (r, s = 1, 2, . . .) with r > s; thus extending the Dobihski 
relations. 

With this in mind we pose the following qüestions: 

(i) What extensions of the conventional Stirling and Bell numbers oceur in the normal 
ordering of [(a"l~) r a s ] n ? 

(ii) Can the generalised Bell numbers B r s {n) so defined be represented by an infinite series 
of the type of Eq.@ - that is, do they satisfy a generalised Dobihski formula ? 

(iii) May one consider the B rjS (n) as the n-th moments of a positive weight function W rjS (x) 
on the positive half-axis, and may this latter be explicitly obtained ? 

In this note we indicate affirmative answers to these qüestions. 
To this end we generalize Eq.@ by defining for r > s: 



k=s 



[(ot)V] n = (at)"^ ]T SV>, k){<J) k a k (10) 
or, alternatively, 

[x r (d/dx) s ] n = x n{r ~ s) S rA n , k ) xk {d/dx) k . (11) 

k=s 

Eqs.([T0|) and < JTT] ) introduce generalized Stirling numbers S r:S (n, k) which imply an extended 
definition of generalized Bell numbers: 

ns 
k=s 
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Note that the Stirling numbers SV,i( n > ^) were studied in [[[Up. Also Bi i(n) = B{n) of Eq.@. 

We have found a representation of the numbers B r ^ s {n) as an infinite series, which is a 
generalization of the Dobihski formula Eq.(]5|). For r = s one obtains: 

n-l 



B rr ( 



n 



'Me) £ 



k=0 



k\ 



(k + r)\ 



n 



1,2 



withS^O) 



1 by convention. For r > s the corresponding formula is: 

\ r—s ' 



B 



n 



[(r 



_ \s(n-l) 



ME 



k=0 



n 



ir(i + S' 



(13) 



(14) 



The formula Eq.(Q) and its extensions Eqs.(fT3|) and ( fÏ4| ) share a common feature, namely, 
the fact that they give rise to a series of integers is by no means evident! 

A general family of sequences arising from Eq.(|ï4|) has the form (p, r = 1, 2, . . .): 



B 



pr+p,pr 



(n) = (1/e) 



n 



(p(n-l)+j)\ 



x r F r (pn + 1, . . . , pn + 1 



x 



p(r - 1); 1 +p, . . . , 1 +p + p(r - 1); 1), (15) 



where r F r is the hypergeometric function. Knowledge of the generalized Stirling numbers in 
Eq.([TO|) solves the normal ordering problem for [(a^) r a s ] n . We are able to give the appropriate 
generating functions for the sequences B r)S (n). It then follows that, at least formally, we can 
furnish the generating functions for S r)S (n, k) as well [JTT]] . Additionally, it turns out that in 
certain circumstances one may obtain explicit expressions for them. We quote two such cases: 

(fe - p - r)!J (r<k<rn) (16) 



and 



Sy y (/FI) /ul 



5*2,1 (n, k) 



k-r 

£ 

p=0 



(k — p)\p\ 



l< k < n) 



(17) 



which are the so-called unsigned Lah numbers [|3|, [K]] . 

For those pairs (r,s) for which we have an explicit expression for S riS (ri) we may 
generalize Eq.(P) to obtain the normal ordered form of e A ^ at ^' aS . For example, the matrix 



element (z|e A(a ^ a \z) leads to the following normally ordered expression: 



3 A(at) r o 



We apply the method to Eq.dTS]) which gave the EGF of Eq.([! 
Eq.@ in the coherent state \z) to get: 

( z \fMr*T*\ z ) =; exp{[(l - \{z*) T ~ x {r - 1))^ 
which evaluates at z = 1 giving 

3 A(ot)'·o 



=: exp{[(l - A(a t ) r - 1 (r - 1))^ - l]a f a} : (18) 

and take the expectation of 



M 2 }: 



(z\e> 



z) z=1 =: exp{(l - X(r - 1))— - 1} : 



(19) 



(20) 



which is precisely the EGF for the numbers B r l (n) |]TQ|] . 

For general r > s the corresponding B r s (n) grow much more rapidly than n! and thus 
may not be obtained via the usual form of EGF One instead defines the EGF in terms 
of B r:S (n) I (n\y where t is an integer chosen to ensure that Y^?=o B r , s (n) / (n\) t+1 has a 
finite radius of convergence. As a result one obtains variants of Eq.<JT9|) involving different 
hypergeometric functions [JTTJ] . 
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The analogue of Eq.@ in the general case 



B r , s {n) = / x n W riS (x)dx 



Jo 



(21) 



leads to weight functions W r>s (x) which may be shown to be positive through use of properties 
of the Mellin transform, and for some of which analytic expressions may be obtained. For 
r = s we obtain discrete measures, giving a "rarefied" form of the Dirac comb, while for 
r > s we obtain continuous measures. We conclude this note with one example of the latter 
kind: 



where l\ (y) is the modified Bessel function of the first kind. 
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